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The goal of this note is to provide most of the technical details involved in the application of 
the Chapman- Enskog method to solve the revised Enskog equation to Navier-Stokes order. Explicit 
expressions for the transport coefficients and the cooling rate are obtained in terms of the coefficient 
of restitution and the solid volume fraction by using a new Sonine approach. This new approach 
consists of replacing, where appropriate in the Chapman-Enskog procedure, the local equilibrium 
distribution (used in the standard first Sonine approximation) by the homogeneous cooling state 
distribution. The calculations are performed in an arbitrary number of dimensions. 



I. FIRST-ORDER DISTRIBUTION FUNCTION 

We consider a granular fluid composed of smooth inelastic hard spheres (d — 3) or disks (d — 2) of mass m and 
diameter a. Collisions are characterized by a (constant) coefficient of normal restitution < a < 1. At a kinetic level, 
all the relevant information on the system is given through the one-particle velocity distribution function /(r, v, t), 
which is assumed to obey the (inelastic) Enskog equation. This equation can be solved by the Chapman-Enskog 
method. To first order in the spatial gradients (Navier-Stokes order), the distribution function /^(r, v, t) is given by 

= A(V) ■ VlnT + B(V) • Vlnn 

+dj (V) i (did, + djUt - %V • v) + V (V) V • U. (1) 



2 V J J d 

The quantities A, 3, and T> are the solutions of the following integral equations [l[ 

1 - ( ;(°)JL.(VA)- 1 -(^A+£A = A, (2) 

i C (0) ^ ■ (VS) + CB = B + C<°> (l + ^ lnx) A, (3) 

\C {0) ^ ■ (VCy) + \C^C %3 + CC %3 = Cy, (4) 

ic (0) ^ • (VP) + \t^V + CV = D, (5) 

where <f> — (Q,d/2 d d)na d is the solid volume fraction and fid = 2iT d ^ 2 /T(d/2) is the area of the d-dimensional unit 
*_h , hard sphere. The linear operator C is given by 

a • 

CX = - (j(°)[/(°),X] + , (6) 

where the collision operator [X, Y] can be recognized as the Boltzmann operator for inelastic collisions multiplied 
by the pair correlation function x(0). The inhomogeneous terms of Eqs. ©-(O are defined by 



1 — "-(0)_ET7, r f(0) 1 



A (V) = -VV V • Vf w - ^Vv/ l,,) + -K V V • ( V/W ) , (7) 
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= _vf(°) - E ( i 
p 



B (V) = -V/ 



^^)vv/^)-(l + ^lnx)/C 



c(0) 



(8) 



(o) 



i f (0) 



(9) 



£> f v • (v/<°>) -Ucu + -/ ] V V • ( V/ 



(0) 



Here, V v = 3/dV, 



and /Ci is the operator 

/C,[X] = * d X (0) I dv 2 / d?8(S ■ g)(a ■ g)a, a - 2 / (0) (vi')^K) + / (0) (vi)^(v 2 ) 



(10) 



(11) 



(12) 



where v'/ = Vi — |(1 + a 1 )((T ■ g)«r and v 2 = v 2 + i(l + a ■ g)S\ In Eq. (|10p . £j/ is defined through the 



expression 



C = C (0) + Cc/V • u, 



(13) 



where is the cooling rate evaluated at zeroth-order. 

With the distribution function / W determined, the goal now is to evaluate the heat and momentum fluxes as well 
as the cooling rate to first order in the spatial gradients. From these constitutive equations one can identify the forms 
of the Navier-Stokes transport coefficients. 



II. PRESSURE TENSOR 



To first order in the spatial gradients, the pressure tensor is given by 



R 



-r, [ diUj + djUi - -8 l3 ■ V • U ) - 7 V • TJS i3 , 



where r\ is the shear viscosity and 7 is the bulk viscosity. Their expressions are given by 

d 



V = Vk 



l + ^ X (l + a) 



d + 2 



7, 



(14) 



(15) 



7 = J B 3 ^-L,^+Wi 1 ~^ 2 - 



Ad 2 



■ma + x(l + a)n«o/ 7 (a). 



(16) 



where n* = na d . The subscript k denotes the contributions from the kinetic part of the pressure tensor. The kinetic 
part 77^. to the shear viscosity is expressed in terms of the solution to the integral equation ^ as 



Vk 



(d- 



T)W^/ dvA ^ (V) ' 



where Dy = m(ViVj — hV Sij). The dimensionless integral L ( is 



J 7 = J_ /dVi / dV 2 /(°)(V 1 )/(°)(V 2 ) 5 , 
n z v J J 



where vo = y/2T/m is the thermal velocity, g = |Vi — V 2 |, and the coefficients B n are defined by 



(17) 



(18) 



(19) 
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The kinetic part r/k follows directly from the integral equation ((4]) by multiplying it by Dij and integrating over the 
velocity. The final result can be written in the form rj k = 770^ , where 



l-!^(l + a)(l-3a)ta 



and 



Vo 



(d + 2)n Y _ d mT 



(20) 



(21) 



is the low density value of the shear viscosity in the elastic limit. In these expressions, Q = (, / v an d I have 
introduced the (reduced) collision frequency 



" ^/dvAi(V)Cy(V)' 

where i/o = nT/r/o. Upon deriving Eq. ([20|) . use has been made of the result 

d 



J dV AjCV)^ 



(0) 



ld-2 



(d-l)nTx«/>(l + Q!)(l-3a). 



(22) 



(23) 



It is worthwhile remarking that, so far the expressions (fT5|) and (|20|) for 77 and (fl6|) for 7 are exact. However, in 
order to get the dependence of both transport coefficients on the coefficient of restitution a and the density <f>, one 
needs to know the explicit form of the reference distribution function and the solution Cu to the integral equation 
(U). The determination of /W to leading order in the Sonine polynomial expansion yields [3j 



f W (V)^f M (V) 1 + - 



m y 2 \ 2 J±o m i/2 



2T 



d + 2 mP d(d + 2) 

~ t~ + 4 



where 



is the local equilibrium distribution function and [3] 

32(l-a)(l-2a 2 ) 



mV 2 /2T 



9 + 24d - a(41 - 8<f) + 30(1 - a)a 2 ' 
Moreover, considering the approximation (|24[) for /(°) and neglecting nonlinear terms in c, the integral 7~ is 



^r(^ 



r(f) 



i-± 

32 



(24) 



(25) 



(26) 



(27) 



With respect to Cij , I'll consider a modified version Q of the Sonine polynomial expansion where the distribution 
f(°'(V) is taking as the weight function in the Sonine expansion instead of the simpler Maxwellian form Jm(V). Thus, 
in the case of the pressure tensor, the leading Sonine approximation to Cjj(V) is 



C«(V) -/ (0) W^TpAi(V). 



(28) 



The collision integral 1/* can be evaluated when one considers the form (|28[) for dj(V). Neglecting nonlinear terms 
in c the result is |4j| 



iiX ( 1 - a + rd ) (1 + a) ( 1 + — c 



4<T 



32 



(29) 
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III. HEAT FLUX 



The constitutive form for the heat flux in the Navier-Stokes approximation is 

q (1) = -kVT - /iVn, 



(30) 



where n is the thermal conductivity, and /i is an additional transport coefficient not present in the elastic case. After 
some algebra, both transport coefficients can be written as 



2 d ~ 2 , . B z mcj d+1 x, s o <,/ 



(31) 



where the dimensionless integral I K is given by 



•>d-2 



1 + 3 TT-S^X( 1 + a ) 



d + 2 



7 - = T3 / dV i / dVa/^CVO/^fVa 



5 - 1 (g.G) 2 + . 9 G 2 + - 9 (g.G) 



and G = i(Vi + V2) is the velocity of center of mass. 
The kinetic parts Kk and ^tfc are defined as 



~ J dvS(V).^V), 



(32) 



(33) 



(34) 



f 1 = 



■i/dvS(V).S(V), 



where 



S,V, = (^-^T)V. 



(35) 



(36) 



The kinetic part of the thermal conductivity is obtained by multiplication of Eq. ([2]) by S(V) and integration over 
the velocity. The result is 



Kk 



1 —{yl-2Q)- x J dVS(V)-A(V) 



where 



dTvo 

^^-^{^^^(l + c)- 1 - I dVS(V)./C Vv(V/(») 



/dvS(V)-£.A(V) 



K v J dvS(V)A(V) ' 

The last term on the right hand side of Eq. (|3T[) can be evaluated more explicitly and the result is 



dVS K 



Vv V/ 



f(0) 



3 n 2 T 2 r c 

= --Sld a d X (l + a) 2 2a-l + -(l + a) 

8 m 12. 



With this result, the kinetic part Kk is given by 



K k = (K ~ 2Co*) l ll + c + 3^-^4>x(l + a) 2 [2a-l + |(l + a) 



where 



K 



d(d + 2) 770 
2{d- 1) m 



(37) 



(38) 



(39) 



(40) 



(41) 
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is the low density value of the thermal conductivity of an elastic gas. 

To evaluate n^, I multiply Eq. by S(V) and integrate over the velocity to get 



dnvQ 



dVS(V)- [(l + ^ln X )C (0) ^+B 



(42) 



where 



/dvS(V)-£(V) 
i/„ /dvS(V) 'B(V)' 



(43) 



The last integral in Eq. flU) is 



d(d + 2) nT 2 / 1 , , 
dVS • B = - A_^_ c _ h + _^ln X 



dvs -k; 



(0) 



(44) 



where the last term can be evaluated in a similar way to Eqs. (|23|) and (|39|) . After some algebra, one gets 

">2 



dvs-x; 



(0) 



= -n d n * Xnr " (l + a) [a(a - 1) + —(10 + 2d - 3a + 3a") 
8 m L 12 



(45) 



With these results, /ik can be written as \j,k = («qT '/ 'n) n* k where 



<f-l 



2 (21/; - 3C *) { Co 4 (1 + 0^ Inx) + -^-c + 3 



2 d ~ 2 {d- 1) 
d(d + 2) 



~ ( 1 + 2^ ln X 



a(a - 1) + —(10 + 2d - 3a + 3a 2 ) 



(46) 



The expressions (pTTj) and (|40p for the thermal conductivity « as well as Eqs. (I3"2l and (|46l) for the coefficient are 
still exact. On the other hand, as in the case of the pressure tensor, to get the explicit forms of k and fi one has to 
determine the collision frequencies v* and v* and the integral I K . As before, the latter integral can be evaluated by 
taking the leading Sonine approximation (|24p for /(°) and neglecting nonlinear terms in c. The result is 



(f) 

With respect to u* and v* one takes now the approximations 



(47) 



A 



1 



d+2l+ dd^cnT 2 



<rn 

K k S(V), B 



m —, s 
MfcS(V), 



where I have introduced the modified Sonine polynomial 



S(V) = S(V) - ^-^TcV. 



(48) 



(49) 



Note that upon writing the approximations (|48[) . I have neglected nonlinear terms in c as well as terms proportional 
to the sixth cumulant of f^ '. The collision integrals i>* and v* can be evaluated with the forms P5|) and the result is 



1 + a 



d-1 3,, „ W1 . 296 + 217d- 3(160 + lld)a ' 
_ + _( d + 8)(l-a) + ^ U 



(50) 



IV. COOLING RATE 



The cooling rate £ is given by 



C - C (0) + Cc/V • u, 



(51) 
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where 

<«._!! + ?(, (52) 

At first order in gradients, the proportionality constant £[/ is a new transport coefficient for granular fluids. This 
coefficient is given by 

(u = (w + Cu, (53) 

where 

Cio = -4^(1 - «Vx = -3^x0(1 - (54) 



Cn = 2^T dr (±ti) - a2 ) J dY ^ J dV 2 .g 3 /W(V 1 )P(V 2 ), (55) 

where the unknown functions 2?(V) are the solutions to the linear integral equation (|5|). 

An approximate solution to the above integral equation (|S|) can be obtained by taking the Sonine approximation 

£>(V) -> e D / M (V)F(V), (56) 

where 

F(V) = (g;)V- ^^ + ^ + 2). (57) 

The coefficient eu is given by 



Substitution of Eq. (|5o| into Eq. (|55|) gives 



3(d+2) ^ 2 ^ 3 



Cu = ^^x(l-^)^l + ^cj, Q e D . (59) 

The coefficient eo is determined by substituting Eq. (|56p into the integral equation ([5]). multiplying by .F(V) and 
integrating over V. The result is 

£ D = (< ~ Iq) 1 -77^— I dVF(V)D(V\ (60) 



7 2 M V d(d+2)ni>b 

where the term £nc has been been neglected in accord with the present approximation. Moreover, the terms propor- 
tional to c coming from v* and Q must be also neglected by consistency. In Eq. (|60p . I have introduced the collision 
frequency 

,_ JdVF(V)C[f M (V)F(V)] 

7 v fdVf M (V)F(V)F(V)- [ ) 

The last integral in Eq. (ISUl) is given by 

i / dVF(V)2>(V) = _^- X (l + a) (i - -)c + -^ TY) /dVF/C, [a«/<°>] , (62) 



d(d + 2) n 

where 



J dVFKt d v J {0) = ^l d n*n X {l + a)\{l-a 2 ){ha-l)- [l5a 3 - 3a 2 + 3(4d + 15)a - (20d + 1)] } 

3 

= g^^d"-*^^*. (63) 
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The last equality of Eq. (|63|) defines the (reduced) quantity ui*. Moreover, the (reduced) collision frequency v* is 
given by 



„ _ 1 + a 
V i ~ 192 



X [30a 3 - 30a 2 + (105 + 24d)a - 56d - 73] 



(64) 



The expression for the coefficient eo can be obtained when one takes into account Eqs. (|52"j) . and (jo"4"|) . The 
result is 



Therefore, the final expression for £[/ can be written as 

'd+2 W3T2) - ( l + a ) (I ~ a ) c 



2 d-2 

~d 



o- j /i 2W 2{d+2) 



64 ^*_rf±2(i_ a 2) 



(65) 



(66) 



V. SUMMARY OF THE MAIN RESULTS 

In this section I provide a summary of the results obtained in this note. First, the reduced hydrostatic pressure p* 



p* = l + 2 d - 2 (l + a )4>x- (67) 
The (reduced) transport coefficients rf — 77 / 770 , 7* = 7A/0, K * — k/kq, and fi* = (n/x/T/co) are given, respectively by 

d . 



»7 =% 



l + |^*x(l + a) 



d + 2 
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(68) 



2 2d+l , 



« — fi 1 



1+ 3^X(l + «) 



(69) 
(70) 



1 + 3^^(1 + a) 



The kinetic parts are 



„* 1 , .* 



v 9 ^° 



-yd-2 



(71) 



(72) 



4 = i^i (*,* - 2C *) 1 |l + c + 3|^^(l + «) 2 [2a-l + |(l + a) 



(73) 



/',. = i " - |Co) {Co4 (1 + Inx) + + 3 2 ^ 2) 1} + a) 

x K + i^lnxJ a(a - 1) + ^(10 + 2d- 3a + 3a 2 ) 

In these expressions, c, Co: ^> aud ^ are functions of a and 0. They are given by 

32(1 - a)(l - 2a 2 ) 



c(a) 



25 + 24d - a(57 - 8d) - 2(1 - a)a 2 ' 



(74) 



(75) 
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c « = ^r (1 - a) n 1+ 32 c 



(76) 



< = ^x(l-a+ld) (l + o) fl + ^c 



32 



(77) 



1 + Q 



d-1 3 , , „ w . 296 + 217d- 3(160 + lld)a 

1 (d+ 8)(1 - a) H — c 

2 16 v A ; 512 



The (reduced) cooling rate £* = C/^o can be written as 



where 



Co 



C=Co+Cc/VV-U, 



(1 - a )y 1 H c , 

4d v M V 32 ' ' 



(78) 



(79) 



(80) 



where 



and 



c^3^Mi-«V d+22(d+2 



- (1 + a) (|- a) 



G4 



3 X 7 8 v x " ; 



1 + a 

192 ' 



[30a 3 - 30a 2 + (105 + 24d)a - 56<i - 73] 



lu* = (1 + a) {(1 -a 2 )(5a- 1) - ^ [l5a 3 - 3a 2 + 3(4d + 15)a - (20d + 1)] } 



(81) 



(82) 



(83) 



Aside from the terms explicitly neglected in Ref. 0], all the above expressions are consistent with those previously 
obtained by J. F. Lutsko Q, except in the coefficient of v* in the denominator of Eq. (fSTj) . Moreover, in the case 
of hard spheres (d = 3), the expressions (|6"8"|) - (l74"]) for the Navier- Stokes transport coefficients and (|81"j) - (f83|) for the 
cooling rate also agree with those obtained in Refs. pj and @, respectively. 
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